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The shape of velocity and temperature profiles near the horizontal conducting plates in 
turbulent Rayleigh-Benard convection are studied numerically and experimentally over 
the Rayleigh number range 10 s < Ra < 3 x 10 11 and the Prandtl number range 0.7 < 
Pr < 5.4. The results show that both the temperature and velocity profiles well agree with 
the classical Prandtl-Blasius laminar boundary-layer profiles, if they are re-sampled in 
the respective dynamical reference frames that fluctuate with the instantaneous thermal 
and velocity boundary-layer thicknesses. 
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1. Introduction 

The turbulent motion in a fluid layer sandwiched by two parallel plates and heated 
from below, i.e. Rayleigh-Benard (RB) convection, has become a fruitful paradigm for 
understanding the physical nature of a wide range of complicated convection problems 
occurring in nature and in engineering problems (Siggia 1994; Ahlers, Lohse & Grossmann 
2009; Lohse & Xia 2010). A key issue in the study of turbulent RB system is to understand 
how heat is transported upwards by turbulent flow across the fluid layer. It is measured 
in terms of the Nusselt number Nu, defined as Nu = J/(kA/H), which depends on the 
turbulent intensity and the fluid properties. These are characterized, respectively, by the 
Rayleigh number Ra and the Prandtl number Pr, namely Ra = agH 3 A/vn and Pr = 
v I k. Here J is the temperature current density across the fluid layer with a height H and 
with an applied temperature difference A, g the gravitational acceleration, and a, and 
k are, respectively, the thermal expansion coefficient, kinematic viscosity, and thermal 
diffusivity of the converting fluid, for which the Oberbeck-Boussinesq approximation is 
considered as valid. As heat transport is controlled by viscosity and thermal diffusion in 
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the immediate vicinity of the solid boundaries, Nu is intimately related to the physics 
of the boundary layers. 

In thermal convective turbulent flow two types of boundary layers (BL) exist near the 
top and bottom plates, both of which are generated and stabilized by the viscous shear 
of the large-scale mean flow: One is the kinematic boundary layer and the other is the 
thermal boundary layer. The two layers are not isolated but are coupled dynamically to 
each other. They both play an essential role in turbulent thermal convection, especially 
for the global heat flux across the fluid layer. Almost all theories proposed to predict the 
relation between Nu and the control parameters Ra and Pr are based on some kind of 
assumptions for the BLs, such a s the stability assumption of the thermal BL from the 
early marg inal stability theory (iMalkusI I1954T). the turb ulent-B L assump t ion from th e 
theories of lShraiman fc Siggial f)l990h and lsiggial f|l994h and of lDubrullej (|200ll . l2002h . 
and the Prandtl-Blasius laminar-BL as s umpt ion of the Grossmann & Lohse (GL) theory 
(jGrossmann fe Lohsell2000l . l200ll 120021 [20o3) . Because of the complicated nature of the 
problem, different theories based on different assumptions for the BL ma y yield the same 
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redactions for the global quan tities, such as the Nu-Ra scaling relation (jCastaing et al 



1989; IShraiman fc Siggialll990f) . Therefore, direct characterization of the BL properties 
is essential for the differences between and the testing of the various theoretical models 
and will also provide insight into the physical nature of turbulent heat transfer in RB 
system. 

In the GL theory, the kinetic energy and thermal dissipation rates have been decom- 
posed into boundary layer and bulk contributions. Scaling wise and in a time averaged 
sense a laminar Prandtl-Blasius boundary layer has been assumed. This theory can suc- 
cessfully describe and predict the N usselt and the Reyn olds number dependences on Ra 
and Pr (see e.g. the recent review in lAhlers et al. 2009). As the Prandtl-Blasius laminar 
BL is a key ingredient of the GL theory, it is important to make direct experimental 
verification of it. We note that also the (experimentally verified) calculation of the mean 
temperature in the bulk i n bot h liquid and gaseous non-Oberbeck-Boussinesq RB flows 
( Ahlers et a71l2006l I2007L 120081 ) is based on the Prandtl-Blasius theory. 

In a recent high-resolution measurement of the properties of the velocity boundary 
layer, Sun, Cheung & Xia (2008) have found that, despite the intermittent emission of 
plumes, the Prandtl-Blasius-type laminar boundary layer description is indeed a good 
approximation, in a time-averaged sense, both in terms of its scaling and its various 
dynamical properties. However, because of the intermittent emissions of thermal plumes 
from the BLs, the detailed dynamics of both kinematic and thermal BLs in turbulent RB 
flow are much more complicated. On the one hand, direct comparison of experimental ve- 
locity (du Puits, Resagk & Thess 2007) and numerical temperature ( Shishkina fe Thessl 
2009l l profiles with theoretical predictions has shown that both the classical Prandtl- 



Blasius laminar BL profile and the empirical turbulent logarithmic profile are not good 
approximations for th e time-averaged velocity and temperature profiles. Furthermore, 
ISugivama et all (|2009h from two-dimensional (2D) and Stevens, Verzicco & Lohse (2010) 
from three-dimensional (3D) numerical simulations found that the deviation of the BL 
profile from the Prandtl-Blasius profile increases from the plate's center towar ds the side- 
walls, due to the rising (falling) plumes near the sidewalls . On the other hand. lQiu fc 
( 19981 ) have found near the sidewall and Sun et al. ( 2008h near the bottom plate that the 
velocity BL obeys the scaling law of the Prandtl-Blasius laminar BL, i.e., its width scales 
as X v /H ~ Re~ 5 , where X v is the kinematic BL thickness, defined as the distance from 
the wall at which the extrapolation of the linear part of the local mean horizontal velocity 
profile u(z) = (u x (z,t)), with z being the vertical distance from the bottom plate and 
(• • •) being the time average at the plate center, meets the horizontal line passing through 
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FlGURE 1. Comparison between the spatial :r- interval and time averaged velocity profiles u(z) 
(triangles), the dynamically rescaled velocity profile u*(zl) (circles - for the notation we refer to 
section 3), and the Prandtl-Blasius velocity profile (solid line) near the bottom plate obtained 
experimentally at Ra — 1.8 x 10 11 and Pr = 5.4 (working fluid water). 



the maximum horizontal velocity [u(z)] max , and Re is the Reynolds number based on 
[u(z)] max . These papers highlight the need to study the nature of the BL profiles, both 
velocity and temperature, in turbulent thermal RB convection. 

Considerable progress on this issue has recently been achieved by I Zhou fc Xia ( 2009h 
who have experimentally studied the velocity BL for water (Pr = 4.3) with particle 
image velocimetry (PIV). They found that, since the dynamics above and below the 
range of the boundary layer is different, a time-average at a fixed height z above the 
plate with respect to the laboratory (or container) frame will sample a mixed dynamics, 
one pertaining to the BL range and the other one pertaining to the bulk, because the 
measurement position will be sometimes inside and sometime outside of the fluctuating 
width of the boundary laye r. To make a clean separation between the two types of 
dynamics, IZhou fc XiaT (|2009l) studied the BL quantities in a time- dependent frame that 
fluctuates with the instantaneous BL thickness itself. Within this dynamical frame, they 
found that the mean velocity profile well agrees with the theoretical Prandtl-Blasius 
laminar BL profile. In figure Q] we show the essence of the results, again for the velocity 
boundary layer but for somewhat larger Pr, now Pr — 5.4. (For details of the experiment 
and the apparatus used, please see Xia, Sun & Zhou 2003; Zhou & Xia 2009). Also here 
th e method of using the time dependent frame works as good as for the Pr — 4.3 case 
of IZhou fc Xial (|2009h . While at the large Ra = 1.8 x 10 11 the time and space averaged 
velocity profile (triangles) already considerably deviates from the Prandtl-Blasius profile 
(solid line), the dynamically rescaled profile (circles) perfectly agrees with the Prandtl- 
Blasius profile. Thus a dynamical algorithm has been established to directly characterize 
the BL properties in turbulent RB systems, which is mathematically well-defined and 
requires no adjustable parameters. 

The questions which immediately arise are: (i) Does this dynamical rescaling method 
also work for the temperature field, giving good agreement with the (Prandtl number 
dependent) Prandtl-Blasius temperature profile? (ii) And does the method also work for 
lower Pr, where the velocity field is more turbulent? Both these questions cannot be 
answered with the current Hong Kong experiments, as PIV only provides the velocity 
field and not the temperature field, and as PIV has not yet been established in gaseous 
RB, i.e., at low Pr number Rayleigh-Benard flows. 
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In the present paper we will answer these two questions with the help of direct numer- 
ical simulations (DNS). To avoid the complications of oscillations and rotations of the 
large scale convection roll plane and as the Pran dtl-Blasius theory is a 2 D theory anyhow 
we will res trict ourselves to th e 2D simulations of lSugivama et al. (2009). Our results will 



show that IZhou &: Xial (|2009l )'s idea of using time-dependent coordinates to disentangle 
the mixed dynamics of BL and bulk works excellently also for the temperature field and 
also for low Pr flow. I.e., if dynamically rescaled, both velocity and temperature BL 
profiles can be brought into excellent agreement with the theoretical Prandtl-Blasius BL 
predictions, for both larger and lower Pr. 



2. DNS of the 2D Oberbeck-Boussinesq equations 

The numerical method has been explained in detail in ISugivama et all ( 20091) . In a 
nutshell, the Oberbeck-Boussinesq equations with no-slip velocity boundary conditions 
at all four walls are solved for a 2D RB cell with a fourth-order finite-difference scheme. 
The aspect ratio is V = D/L = 1.0, the Rayleigh numb er Ra = 10 8 — 10 9 , an d the 
Prandtl number either Pr — 4.3 (water) or Pr — 0.7 (gas). Sugivama et al. ( 20091 ) have 
provided a detailed code validation. 

As the governing equations are strictly Oberbeck-Boussinesq, there exists a top-bottom 
symmetry. We therefore discuss only the velocity and temperature profiles near the bot- 
tom plate. For the temperature profiles, we introduce the non-dimensional temperature 
0(z, t), defined as 

e M = -^M, (2-i) 

where 9 bot is the temperature of the bottom plate. In this definition, <d(H) = 2 and 
O(0) = are the temperatures for the top and bottom plates, respectively, and Q(H/2) = 
1 is the mean bulk temperature. 



3. Dynamical BL rescaling 

The idea of the IZhou fc Xial (120091) method is to construct a dynamical frame that 
fluctuates with the local instantaneous BL thickness. To do this, first the instantaneous 
ki nematic and t herma l BL thicknesses are determined using the algorithm introduced 
by IZhou &: Xia ( 20091 ). To reduce data scatter, the horizontal velocity and temperature 
profiles at each discrete time t, u(z, t) and Q(z, t), are obtained by averaging the velocity 
and temperature fields along the x-direction (horizontal) over the range 0.475 < x/D < 
0.525. Figures [2ja) and (b) show examples of u(z,t) and <d(z,t) versus the normalized 
height z/H, respectively, of the DNS data obtained at Ra — 10 9 and Pr = 0.7. Both 
u(z, t) and 0(z, t) rise very quickly from to either the instantaneous maximum velocity 
or to the bulk temperature within very thin layers above the bottom plate. While after 
reaching its maximum value, u(z,t) slowly decreases in the bulk region of the closed 
convection cell, Q{z,t) reaches and stays nearly constant at the bulk temperature = 1. 
To see the velocity and the temperature in the vicinity of plates more resolved, we plot 
the enlarged near-plate parts of the u(z,t) and Q(z,t) profiles in figures [2] (c) and (d). 
One sees that both profiles enjoy a linear portion near the plate. The instantaneous 
velocity BL thickness 8 v (t) is then defined as the distance from the plate at which the 
extrapolation of the linear part of the velocity profile meets the horizontal line passing 
through the instantaneous maximum horizontal velocity, and the instantaneous thermal 
BL thickness Sth(t) is obtained as the distance from the plate at which the extrapolation 



Prandtl-Blasius temperature and velocity BL profiles in turbulent RBC 



o.io 



o.oo - 



-o.io 




0.0 0.2 0.4 0.6 


0.8 1.0 


z/H 




'/ 1 1 
/ 

/ 


1 




* oeeoooooooD °Obo 


i 
i 

! 


(c) 


' , 1 


1 



0.00 



z/H 



2.0 



1.5 



1.0 - 



0.5 - 



0.0 



n — i — i — i — i — i — i — i — i — i — r 

(b) 




J i I i I i I i I i L 

0.0 0.2 0.4 0.6 0.8 1.0 

Z/H 



0. 



^ 0.6 I— 
.4 I- / 



0.2 



/ 



1.0 l 

/ 

// 




°Q*?OQfiPQQOOOOOOOOOOOO< 



J 



(d) 



0.00 



z/H 



Figure 2. Examples of (a) an instantaneous horizontal velocity profile u(z,t) and (b) a normal- 
ized instantaneous temperature profile G(z, i), averaged over 0.475 < x/D < 0.525. The DNS 
data are obtained at Ra = fO 9 and Pr — 0.7. (c) and (d) show enlarged portions of the velocity 
and temperature profiles near the bottom plate, respectively. The two tilted dashed lines are 
linear fits to the linear parts of the velocity and temperature profiles near the plate and the 
two horizontal dashed lines mark the instantaneous maximum horizontal velocity and the bulk 
temperature O = f , respectively. The distances of there crossing points from the plate define the 
instantaneous BL thicknesses y,(t). The instantaneous profiles are not top-down symmetric, 
the time averaged ones are. Within our present statistical error our data are consistent with 
zero thermal gradient in the bulk. 



of the linear part of the temperature profile crosses the horizontal line passing through 
the bulk temperature. The arrows in figures [2jc) and (d) illustrate how to determine 
5 v (t) and S t h(t) as the crossing point distances. 

With these measured 8 v (t) and Sth(t), we can now construct the local dynamical BL 
frames at the plate's center. The time-dependent rescaled distances z*(t) and z% h (t) from 
the plate in terms of 5 v (t) and 5th{t), respectively, are defined as 

z*(t) =z/S v (t) and z* th (t)=z/5 th (t). (3.1) 

The dynamically time averaged mean velocity and temperature profiles u*(z*) and Q*(z^ h ) 
in the dynamical BL frames are then obtained by averaging over all values of u(z, t) and 
Q{z,t) that were measured at different discrete times t but at the same relative positions 
z* and z* h , respectively, i.e., 

u *( z v) = {u(z,t)\z = z*S v (t)) and Q*(z* h ) = (Q(z,t)\z = zl h 5 t h(t)). (3.2) 

We first discuss our results from the simulation performed at Pr — 4.3, the Prandtl 
number corresponding to water at 40 °C. Figure[3Ja) shows the u*(z*) profile (circles), 
normalized by its maximum value [u*(z*)] max , obtained at Ra = 10 s . For comparison, 
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Figure 3. Comparison among (a) velocity profiles: dynamical frame based u*(z*) (circles), 
laboratory frame based u(z) (triangles), and the Prandtl-Blasius velocity profile (solid line), 
and (b) the corresponding temperature profiles: Q*(z^ h ) (circles), Q(z) (triangles), and the 
Prandtl-Blasius temperature profile (solid line) near the bottom plate. All results obtained 
numerically at Ra = 10 8 and Pr — 4.3. The inset of (b) shows enlarged portions of the profiles 
around the thermal boundary layers' mergers to the bulk. 



we also plot in the figure the time-averaged horizontal velocity profile u(z) (= (u(z,t)}) 
(triangles), obtained from the same simulation. The solid line represents the Prandtl- 
Blasius velo city BL profile, the initial slope of which is matched to that of the measured 
profiles (cf. lAhlers et all [20 0(f). For the range z* < 2 the u*(z*) profile obtained in the 
dynamical frame agrees well with the Prandtl-Blasius profile, while the time-averaged 
u{z) profile obtained in the laboratory frame obviously is much lower than the Prandtl- 
Blasius profile in the region around a few kinematic BL widths. - Note that for z* > 2 
the u*(z*) profile deviates gradually from the Prandtl-Blasius profile because u*(z*) 
decreases in the bulk region of the closed convection system down to in the center 
and then changes sign. The Prandtl-Blasius profile, instead, describes the situation of an 
asymptotically constant, nonzero flow velo city. - These DNS results are similar to those 
found experimentally in a rectangular cell (jZhou fc Xlall2009h . 

Figure [3Jb) shows a direct comparison among the temperature profiles obtained from 
the same simulation: the dynamical frame based Q*(z* h ) (circles), the laboratory frame 
time-averaged temperature profile Q(z) (= (Q(z, t))) (triangles), and the Prandtl-Blasius 
temperature profile. At first glance both the Q*(z^ h ) and &(z) profiles are consistent with 
the Prandtl-Blasius thermal profile. However, looking more carefully at the region around 
the thermal BL to bulk merger (the inset of figure^b)), one notes that the Q*(z^ h ) profile 
obtained in the dynamical frame is significantly closer to the Prandtl-Blasius profile 
than the time-averaged 8(z) profile obta i ned i n the laboratory frame, indicating that 
the dynamical frame idea of I Zhou fc Xial (2009) works also for the thermal BL. Taken 
together, figures a) and (b) illustrate that both the kinematic and the thermal BLs in 
turbulent RB convection are of Prandtl-Blasius type, w hich is a key assumption of the 
G L theory dGrossmann fc Lohsel[2000l l200lL I2OO2I |2004> . and the dynamical frame idea 
of I Zhou &: Xial (|2009l ) can achieve a clean separation for both temperature and velocity 
fields between their BL and bulk dynamics. 

We next turn to the simulation performed at Pr — 0.7, a Prandtl number typical for 
gases, which is relevant in all atmospheric processes and many technical applications. 
Figures Ufa) and (b) show direct comparison between the temperature and velocity pro- 
files, respectively, at Ra = 10 9 . Again, around the BL-bulk merger range the laboratory 
frame time-averaged profiles are found to be obviously lower than the Prandtl-Blasius 
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Figure 4. Comparison between (a) velocity profiles: dynamical u*(z*) (circles), laboratory u(z) 
(triangles), and the Prandtl-Blasius laminar velocity profile (solid line), and (b) temperature 
profiles: dynamical Q*(Zf h ) (circles), laboratory 0(z) (triangles), and the Prandtl-Blasius lami- 
nar temperature profile (solid line) near the bottom plate, all obtained numerically at Ra — 10 9 
and Pr = 0.7, representative for gases. 



profile. This once more indicates that the time-averaged BL quantities obtained in the 
laboratory frame are contaminated by the mixed dynamics inside and outside the fluctu- 
ating BLs. On the other hand, within the dynamical frame, both u*(z*) and Q*(z^ h ) are 
found to agree pretty well with the Prandtl-Blasius laminar BL profiles, indicating that 
the dynamical frame idea works also for the turbulent RB system with working fluids 
whose Prandtl numbers are of the same order as those for gases. 



4. Shape factors of the velocity and temperature profiles 

Let us now quantitatively compare the differences between the Prandtl-Blasius profile 
and the profiles obtained from both simulations and experiments for various Ra and 
various Pr. The shapes of the velocity and temperature (thermal) profiles, labeled by 
i '. = v or i = th, can be chara cterized quantitatively by their shape factors Hi, defined as 



(jSchlichting fc Gerstenll2004D 



Hi = -±- , i = v,th. (4.1) 

Xf and A™ denote, respectively, the displacement and the momentum thicknesses of the 
profile, namely, 

1 [Y(z)} 111 (IX max 

Here Y(z) — u(z) is the velocity profile if i = v and Y(z) = O(z) the thermal profile if 
i = th. The deviation of these profiles from the Prandtl-Blasius profile is then measured 

by 

5 Hi = Hi — Hf B , (4.3) 

where Hf B is the shape factor for the respective Prandtl-Blasius laminar BL profile. If 
a given profile exactly matches the Prandtl-Blasius profile, SHi is zero. Note that the 
Prandtl-Blasius velocity profile shape factor H BB — 2.59 is independent of Pr, while the 
thermal Prandtl-Blasius BL profile shape factor Hf h B varies with Pr. 

Figure EJa) shows the shape factors Hi(Pr) of the thermal and the velocity Prandtl- 
Blasius BL profiles as functions of Pr and figure EJb) shows the corresponding thermal 
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Figure 5. (a) The shape factors for the thermal (solid line) and velocity (dashed line) 
Prandtl-Blasius BL profiles as function of Pr. The asymptotic value H th (Pr S> 1) is 2.61676... 
and Hv B = 2.59. Both Prandtl-Blasius BL profiles for the velocity and for the temperature 
coincide for Pr = 1. (b) The thermal Prandtl-Blasius BL profiles for three (four) Pr numbers 
and the reference linear and exponential profiles, see (|4.4jl ; in the figure's resolution the thermal 
profile for Pr = 100 is indistinguishable from that of Pr = 1. Note that the shape factor of the 
thermal Prandtl-Blasius BL profile decreases with decreasing Pr due to the slower approach to 
its asymptotic level 1. 



profiles as functions of z* th for three different Pr. Note that the Prandtl-Blasius velocity 
BL profile is identical to the thermal one for Pr = 1. The two figures show that the ther- 
mal shape factor H^ B decreases with decreasing Pr. We attribute this to the decrease 
of the temperature profiles in the BL range and the corresponding increase of the tails 
for lower Pr. Thus we expect that the slower approach to the asymptotic height 1 of 
the thermal profiles in the laboratory frame in figures |3] and U should lead to a nega- 
tive deviation of their H t h's from the respective Prandtl-Blasius values, cf. figure [6l In 
contrast, a positive SHi is obtained if the profile runs to its asymptotic level faster than 
the Prandtl-Blasius profile. To see this more clearly, we have plotted in figure EJb) also 
two extreme cases, the linear and the exponential profiles. Using (|4.2[) one calculates the 
shape factor 3 for the linear profile Q*(z^ h ) — min(l, z^ h ) and the shape factor 2 for the 
exponential one Q*(z^ h ) = 1 — cxp(— z* h ). The iJ-decreasing effect by lowering the profile 
can also be demonstrated by analyzing some profiles analytically. Using a combination 
of exponential profiles, 

8*04) = 0.5(1 - exp(-(l - n){z* th )) + 0.5(1 - exp(-(l + n)z* th )), (4.4) 

with ^ n < 1 one can evaluate, using (|4.2p . that the shape factor for small n is 

H(n) » H(n = 0) - n 2 = 2 - n 2 . (4.5) 

As is shown in figure [5][b) the profile for n > is below the profile for n = 0. This 
analytical example again reflects what we found as the characteristic difference between 
the laboratory frame profiles as compared to the dynamical frame profiles. 

Figure EJa) shows the velocity shape factor deviations SH V (open symbols) and 5H* 
(solid symbols) as obtained from simulations at Pr — 0.7 (circles) and Pr = 4.3 (trian- 
gles) as well as from experiments at Pr = 5.4 (squares). Here, SH V is calculated with 
the time averaged profile u(z) in the laboratory frame, while 5H* is calculated with the 
dynamical, time-dependent frame profile u*(z*). The laboratoty frame based deviations 
turn out to be definitely smaller than zero. In contrast, the shape factor deviations SH* 
for the dynamical frame profiles obviously are much closer to zero. A similar result is 
found for the thermal BLs: Figure EJb) shows SH t h (open symbols) and SHf h (solid 
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symbols), versus Ra, for the same Pr number simulations. Again 8H^ h is nearly zero, 
whereas 8H t h is significantly off. Thus these quantitative deviation measures again indi- 
cate that the algorithm using the dynamical coordinates can effectively disentangle the 
mixed dynamics inside and outside the fluctuating BLs. 



5. Conclusions 

In summary, we have studied the velocity and temperature BL profiles in turbulent 
RB convection both numerically and experimentally. We extended previous results to 
different Prandtl numbers and in particular to thermal BLs. The results show that both 
the velocity and the temperature BLs (at least in the plates' center region) are of laminar 
Prandtl-Blasius type in the co-moving dynamical frame in turbulent thermal convection 
for the parameter ranges studied. However, the fluctuations of the BL widths, induced 
by the fluctuations of the large-scale mean flow and the emissions of thermal plumes, 
cause measuring probes at fixed heights above the plate to sample a mixed dynamics, 
one pertaining to the BL range and the other one pertaining to the bulk. This is the 
reason why the time-averaged velocity and temperature profiles measured in previous 
work in fixed laboratory (RB cell) frames deviate from the Prandtl-Blasius profiles. To 
disentangle that mixed dynamics, we constructed a dynamical BL frame that fluctuates 
with the instantaneous BL thicknesses. Within this dynamical frame, both velocity and 
temperature profiles are very well consistent with the classical Prandtl-Blasius laminar 
BL profiles, both for lower and larger Pr (from 0.7 to 5.4). We have thus validated the 
idea and algorithm of using dynamical coordinates over a range of Pr and Ra for both 
kinematic and thermal BLs and have shown that the Prandtl-Blasius laminar BL profile 
is a valid description for the BLs of both velocity and temperature in turbulent thermal 
convection. Laminar Prandtl-Blasius BL theory in turbulent RB thermal convection has 
thus turned out to indeed be valid not only scaling wise, but also in the time average as 
seen from the dynamical frame, co-moving with the local, instantaneous BL widths. 
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